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ABSTRACT 

We have searched for IR emission from the intracluster dust (ICD) in the galaxy cluster A2029. 
Weak signals of enhanced extended emission in the cluster are detected at 24 and 70 /im. However, 
the signals are indistinguishable from the foreground fluctuations. The 24 versus 70 jxm color map 
does not discriminate the dust emission in the cluster from the cirrus emission. After excluding the 
contamination from the point sources, we obtain upper limits for the extended ICD emission in A2029, 
5 x 10 3 Jy sr -1 at 24 /im and 5 x 10 4 Jy sr -1 at 70 /zm. The upper limits are generally consistent with 
the expectation from theoretical calculations and support a dust deficiency in the cluster compared to 
the ISM in our galaxy. Our results suggest that even with the much improved sensitivity of current 
IR telescopes, a clear detection of the IR emission from ICD may be difficult due to cirrus noise. 
Subject headings: galaxies: clusters: individual: (Abell 2029) — galaxies: clusters: general — inter- 
galactic medium — infrared: general 



1. INTRODUCTION 

An important constituent of the X-ray-emitting gas 
in massive clusters is thought to be material lost from 
cluster galaxies. Such gas would be well processed and 
contain a significant component of heavy elements, indi- 
cated by the absorption lines in the X-ray spectra. This 
galactic-originated gas would also introduce dust grains 
into the intracluster (IC) gas, which can survive sput- 
tering by ambient gas for a typical time scale of ~ 10 8 
yr. However, in contrast to the heavy elements, solid 
evidence of the presence of dust grains in IC gas is still 
lacking. The indirect measurement of ICD, suggested 
by enhanced visual extinction towards clusters deduced 
from high-r edshift objects, is still very controversial (e.g., 
[MaLlfl995h . 

A direct measurement of ICD can be provided by its IR 
emission when they are col li sional ly heated by hot X-ray 
emitting gas. iDwek et al.1 (jl.990) modeled this process 
and predict flux levels of about 0.2 MJy sr -1 at 100 /im 
in the Coma cluster. Their calculations, and the upper 
limit constraints provided by IRAS suggest a deficiency 
of dust in the central region by a factor of 100 compared 
with the interstella r dust-to-gas ratio. Using ISO data, 
IStickel et al.1 (|1998( ) report the detection of a flux excess 
of ~ 0.1 MJy sr -1 at 120 /im in the central region of the 
Coma cluster and a t tribut e this excess to the ICD emis- 
sion. IStickel et all (|2002f) confirmed this excess in the 
I120 /im/^iso /itm color profile across the Coma core re- 
gion, but they did not find the same signa l in the other 5 
clusters they studied. lOuillen et alj (|1999f ) estimated the 
flux in Coma that might arise from faint infrared-excess 
galaxies and concluded that it could be a signi ficant frac- 
tion o f the apparent extended component. iBai et aTl 
( 2006) have studied the Coma cluster and found the pro- 
jected number density of infrared-emitting galaxies is the 
highest in the core region. Such an excess of IR galaxies 
could cause an excess of IR emission in the cluster core 
region, raising more questions about the origin of the 
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excess emission reported by IStickel et al.l (fl99 8. 2002). 

The observational situation is therefore muddled. The 
difficulties arc mostly associated with the sensitivity and 
detector stability in previous missions, but also with their 
limited angular resolution which has made it difficult to 
distinguish the composite of the emission from individual 
cluster galaxies from that of ICD. The low resolution also 
causes confusion limitations due to extragalactic sources. 
But with the much-improved sensitivity and spatial res- 
olution of current and future IR telescopes, e.g., Spitzer, 
AKARI and Herschel, we may marginally reach the de- 
tection thresholds predicted by theoretical calculations 
(|Yamada fc Kitavamall2005h . 

In this paper, we present an attempt to measure the 
extended IR emission from ICD in the galaxy cluster 
A2029 at z = 0.077. It is one of the brightest X-ray clus- 
ters nearby, with a X-ray luminosity Lx(2 — lOkeV) = 
1.1 x I0 45 hjy ergs s" 1 (jDavid et al.lll993l ) and an intra- 
cluster gasmetajlicity about a third of the solar metal- 
licity () White) 12000). The cluster is in a region of rela- 
tively low IR cirrus, which makes it a perfect target for 
searching for the ICD emission. The size of its X-ray 
emitting gas (~ 2'), which limits the size of the ICD 
emissio n, is well suited to s ensitive measurements with 
Spitzer. [Hansen et al.l (|2000f l attempted to assess the IR 
emission from ICD in this cluster with ISO data, but 
with no success. Using MIPS observations, we can dis- 
entangle the IR emission of the galaxies from that of the 
ICD in A2029 and are able to set upper limits to the 
ICD emission at 24 and 70 /im. Throughout the pa- 
per, we assume a ACDM cosmology with parameter set 
(A,n ,Ao) = (0.7,0.3,0.7). 

2. OBSERVATIONS AND DATA REDUCTION 

The cluster was first observed in Feb. 2004 with MIPS 
in a medium rate scan mode. The 24 /im, 70 /im and 
160 /zm data were acquired simultaneously for the central 
40' x 60' region. The integration time was about 80, 40, 
and 8 seconds at 24, 70 and 160 /im. To further confirm 
the weak extended emission detected in the cluster core 
region, a second observation was obtained in Jan, 2005 
to tripled the observation time. 
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The data were proces sed with the MIPS D ata Analysis 
Tool (DAT version 2.9: iGordon dTall 12005ft . At 24 /im, 
extra steps were carried out to remove scattered light, 
which appears as a low-level background modulation syn- 
chronous with scan mirror position. We generated aver- 
age backgrounds for each scan mirror position and sub- 
tracted them from the images. The final mosaic com- 
bining old and new observations has an exposure time of 
~ 340 s pixel -1 on average. For 70 /im, we only use the 
data from the second observation because the first obser- 
vation was taken before the bias was changed on the 70 
/im array and has quite noticeable stim latents. The to- 
tal exposure time for the 70 /im image is ~ 120 s pixel -1 
and is ~ 30 s pixel -1 for the 160 /im image. At the 70 
and 160 /im bands, we took extra processing steps to 
improve the images. The pixel-dependent backgrounds 
were subtracted from the images, which were obtained 
by fitting the images outside of the central 6' region of 
the cluster. This process tends to de-emphasize the ex- 
tended emission outside the source region, but it keeps 
the signal we are interested in. The mosaics have point- 
spread functions (PSF) with FWHMs of ~ 6", 18", 40" 
at 24, 70 and 160 /im. 

To avoid contamination from individual galaxies, we 
extract all the IR sources from the images using Sex- 
tractor. For the 24 /jm image, we exclude all the pixels 
contributing to IR sources down to the 1 a detection limit 
(~ 45 (J,Jy). After the source extraction, 74% of the pix- 
els remain in the 24 /im image. For the 70 and 160 /im 
data, due to the lower uniformity of the mosaic image, 
we only exclude the IR sources down to the 2 a (~ 7 mJy 
and ~ 50 mJy, respectively) sensitivity level, after which, 
76% and 85% of the pixels remain at each band. If the 
IR source density in the cluster core region were greatly 
enhanced by the cluster galaxies, this source extraction 
process might cause a bias when we try to measure the 
extended emission from the remaining pixels. However, 
this is not the case. The fractions of the remaining pix- 
els in the central r < 2' region of the cluster are 70%, 
80% and 90% at 24, 70 and 160 /im, respectively, which 
do not differ significantly from the average values for the 
whole images. 

3. RESULTS 

After excluding the IR sources, we smooth the images 
using a median filter with a box size of ~ 2.'5. The filter 
matches roughly the size of the X -ray emitting regio n (to 
~ 10% of the central gas density. [Lewis et aLll2003ft . All 
IR images show significant brightness fluctuations across 
the map, and the patterns of the fluctuations are roughly 
consistent at the three wavelengths. The fluctuations are 
probably due to the foreground cirrus emission, but they 
could also be caused by the fluctuation of the number 
density of unresolved faint sources, especially at 70 /im, 
where the observations have not reached the conf usion 
limit due to extragalactic sources (|Dole et al.ll2004f) and 
the cirrus emission is still relatively weak. 

In Fig. [I] we show the smoothed images with the X- 
ray contours overlaid. The X-ray data has an expo- 
sure time of 79 ks and were retrieved from the Chan- 
dra Data Archive (ObsID 4977). The 24 /im smoothed 
image shows a weak enhancement corresponding to the 
X-ray emitting region, but it is indistinguishable from 
the foreground cirrus structure and its morphology has 



no apparent similarity with the X-ray image. At 70 /im, 
the enhancement is even harder to recognize due to a 
large cloud of faint emission crossing the cluster. The 
core region is also slightly fainter than the outer region 
of the X-ray emitting region. The enhancement, if any, is 
likely just to be the extended emission from the Galactic 
cirrus. At 160 /im, the X-ray emitting region is totally 
indistinguishable from the b ackground. 

Using ISO observations, iStickel et all (|1998l 12002ft 
found a peak in the /120 /im/^iso /im surface brightness 
ratio towards the Coma cluster center and they use this 
color difference to discriminate the ICD emission from 
the foreground cirrus. To see if a similar enhancement 
of the brightness ratio occurs in A2029, we derived a 
color map by dividing the 24 /im image by the 70 /im 
image. Since the backgrounds have already been sub- 
tracted from both images and we are mostly interested 
in comparing the color of the cluster emission with that 
of the interstellar cirrus, we add the average brightness 
of the interstellar cirrus estimated by SPOT 2 before we 
calculate the brightness ratio. In Fig. [21 the X-ray emis- 
sion region shows a warmer color than the faint cloud 
shown in the 70 /im map. However, again, this enhance- 
ment is not significant compared to the fluctuations in 
the color map, and there are some other patches in the 
image showing similar color to the X-ray emission region. 
We do not show the /70 /jm/Iwo /im color map here be- 
cause the color is dominated by the fluctuations at 70 /im 
due to the much larger cirrus emission at 160 /im than 
at 70 /im. 

Although the extended IR emission of the cluster is in- 
distinguishable from the foreground cirrus in the images, 
we can obtain an upper limit to the ICD emission. In 
Fig. [31 we plot the average surface brightness as a func- 
tion of distance from the cluster center, derived from 
the point sources-removed, non-smoothed image. The 
error bars are the standard error of the mean for each 
data point. The 24 /im radial profile shows a clear peak 
at radius < 1', coincident with the peak of the X-ray 
brightness profile. The enhancement within r < 1' rela- 
tive to the average of the whole image sets a 2a upper 
limit of 5 x 10 3 Jy sr _1 for the extended emission from 
the ICD. However, the apparent significance of the peak 
in the radial profile is mostly due to the fact that at large 
radii, the surface brightness is the average of pixels in a 
larger area, in which the small scale flux density fluc- 
tuations have already been canceled out. To show the 
fluctuation of the brightness at a smaller scale, we plot 
the histogram of the pixels in the IR sources masked-out, 
smoothed 24 /im image. The smoothing uses a median 
box of size 80"; it retains the fluctuations at spatial scales 
larger than the box size. The surface brightness at the 
center of the cluster, although also showing an enhance- 
ment (~ 4 x 10 3 Jy sr -1 ), is only la higher than the 
average of the smoothed image. 

At 70 /im, the enhancement in the cluster is less ob- 
vious due to the foreground cloud it lies in. The peak 
of the emission is not at the center of the cluster, but 
about 2' away. Therefore we estimate the enhancement 
in the central r < 2' region, and obtain a 2a upper limit 
of 1 x 10 5 Jy sr -1 . However, if we exclude the emis- 

2 SPOT, the Spitzer tool for planning observations and submit- 
ting proposals, http://ssc.spitzer.caltech.edu/propkit/spot/ 
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sion in the cluster coming from this foreground cloud 
and compare the enhancement of the cluster relative to 
the local average brightness in the r < 15' region (ap- 
proximately the scale of the cloud), we have a smaller 
upper limit, 5 x 10 4 Jy sr _1 . The cloud also causes the 
histogram of the smoothed image (with a median box 
of size 160") to peak slightly off the mean. Although 
the surface brightness at the cluster center shows an en- 
hancement at about 0.9<r relative to the average, the 
enhancement is less significant compared to the peak 
of the distribution. We can also compare our results 
with the average brightness fluctuation o f cirrus emission 
given by iMiville-Deschenes et al.1 ()2007l ). which is about 
2 x 10 4 Jy sr _1 at this scale for a low brightness region. 
The upper limit we measured is only a few times larger 
than the average cirrus fluctuation, which again proves 
the insignificance of the enhancement compared to cirrus 
noise. At 160 /im, the cluster region shows a lower than 
average surface density compared to the whole image. 
The 2cr error of this non-detection is 6 x 10 4 Jy sr~ 4 . 

4. DISCUSSION AND CONCLUSIONS 

iDwek et al.l (|1990t ) calculated the infrared emission of 
the ICD collisionally heated by X-ray-emitting hot gas 
in the Coma cluster, which has a similar temperature 
and me tallicity, but a higher density than the hot gas in 
A2029 (|Whitell2000T ). Therefore the expected dust emis- 
sion in the Coma cluster should be larger than that in 
A2029 and can provide an upper limit to the dust emis- 
sion in A2029. In Coma, at a distance of about 100 
kpc from the cluster center, the expected dust emission 
is about 5 x 10 2 Jy sr~ 4 at 25 /im and 3 x 10 4 Jy sr _1 
at 60 /im. The upper limit we obtained at 24 /an, 
5 x 10 3 Jy sr- 1 within -90 kpc of A2029, is 10 times 
larger than the expected intensity. At 70 /im, the up- 
per limit excluding emission from the foreground cloud, 
5 x 10 4 Jy sr _1 , which does not change much from r < 90 
kpc to r < 180 kp c, is very clo s e to t he expected inten- 
sity. According to IDwek et al.1 (|1990l ). their calculation 
suggests the dust in the central region of a cluster is defi- 
cient by a factor of 100 relative to the standard value for 
the ISM. The agreement between our upper limit mea- 
surements and their calculation confirms that the dust, 
if any, should be at leas t defic ient by that amount. 

lYamada fc Kitavamal (|2005| ) performed a comprehen- 
sive study of the emission of ICD in a sample of clusters 
with z — 0.01 — 0.8. They present the expected 70 /im 
intensities at the projected distance of 20 kpc from the 
center of the clusters (Fig. 7 of their paper) . In order to 
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compare their results with the upper limits we obtained 
for A2029, we scale the intensities at 20 kpc with a factor 
of 0.28 to deduce the average intensities in the central 175 
kpc region, the same region in which we measured the up- 
per limits for A2029. The factor is obtained by assuming 
the clusters have a similar intensity radial profile as illus- 
trated by the Perseus Cluster in their paper. The results 
are shown in Fig. [4j The two upper limits we obtained 
are plotted as the downward arrows. They are very close 
to the predicted intensities of the clusters at similar red- 
shift (the large circles and triangles). In particular, there 
is one cluster, A3112 (shown as the filled triangle), with 
very similar properties to A2029. They both are regular, 
X-ray bright clusters with redshift z ~ 0.07 and at sim- 
ila r Galac t ic lat itudes. Ac cording to the X-ray analysis 
bv I White! (|2000l ). on which lYamada fc Kitavamal (|2005l ) 
base their calculation, A2029 has very similar metallicity, 
gas temperature and density compared to A3112, which 
would indicate an IR ICD emission at the same level as 
that in A3112. The expected intensity in A3112 is very 
close to the conservative upper limit we measured but is 
higher than the upper limit that excludes the emission 
from the foreground cirrus. The discrepancy suggests 
that t he theoretical calculation bv lYamada fc Kitavamal 
( 2005) may overestimate the observed ICD emission. 

Another way to detect the weak IR e mission from ICD 
is by stacking many galaxy clusters. iMontier fc Giardl 
(2005) stacked a large sample (11,507) of IRAS cluster 
maps together and found a statistical detection of 2.1 ± 
0.7 xlO 3 Jy sr^ 1 at 25 /im and 1.6±0.7xl0 4 Jy sr" 1 at 60 
/im, both several times smaller than the upper limits we 
obtained in A2029. This is consistent with the fact that 
A2029 is one of the brightest X-ray clusters and we would 
expect the ICD emission in this cluster to be higher than 
the average value of a large sample of clusters. 

The faint average value from stacking plus our failure 
to detect significant signals in a cluster expected to be 
relatively bright together indicate that the IR emission 
of ICD gas must be generally very faint. Even with the 
much improved sensitivities of the current generation of 
IR telescopes, this faint ICD emission will be very hard 
to be distinguished from cirrus fluctuations. 
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Fig. 2. — The surface brightness map of 724 /zmi with X-ray contours overlaid. The size and direction of the image, as well as 

the scale of the contours, are the same as in Fig. [2] 
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Fig. 3. — Top: the radial profiles of brightness at 24 and 70 /mi as a function of distance from the center of the cluster. The error bars 
are the errors of the mean in each region. The dashed lines are the average values for each image. The X-ray brightness profiles are plotted 
as the dotted lines. Bottom: the histograms of the smoothed images. The 24 /mi image is smoothed with a median box of size 80" and 
the 70 fim image is smoothed with a median box of size 160". The vertical lines are the values at the center of the cluster. 
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Fig. 4. — Expected 70 /an intensities of ICD in the central 175 kpc of a sample of clusters from Yamada & Kitayama (2005) compared 
with the upper limits measured in A2029 in the same region. The triangles are for the five clusters with highest 70 (im intensity. The filled 
triangle is A3112, which we expect to have a similar level of emission as A2029. The large circles are clusters at rcdshift 0.01 < z < 0.1, 
medium circles are at 0.1 < z < 0.3 and small circles are at 0.3 < z < 0.8. The top downward arrow is the 2rr upper limit for A2029 
measured relative to the whole 70 /mi image; the lower downward arrow is the upper limit obtained by subtracting the local foreground 
emission (r < 15'). 



